Glass Phenomenology from the Connection to Spin Glasses 
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Typical features of glass phenomenology such as the Vogel-Fulcher law, the Kauzmann paradox 
and the Adam-Gibbs relationship are shown to follow from the recently discovered mapping of glasses 
to Ising spin glasses in a magnetic field. There seems to be sufficient universality near the glass 
transition temperature T g such that study of the spin glass system can provide semi-quantitative 
results for supercooled liquids. 
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I. INTRODUCTION 

Under fast enough cooling or densification, materials 
which are as diverse as molecular and polymeric liquids, 
colloidal suspensions, granular assemblies and molten 
mixtures of metallic atoms, may form glasses [lj. These 
are amorphous states that may be characterized mechan- 
ically as a solid, but lack long-range crystalline order. 
Despite all the work devoted to the subject, the mecha- 
nisms responsible for the vitrification processes are not 
well understood, and the transition to the glassy state 
remains one of the most controversial problems in con- 
densed matter physics. 

In a recent paper 0], using an effective potential 
method, a replica formalism has been set up to describe 
supercooled liquids. This approach shows that the ther- 
modynamics of these systems near their glass transition 
temperature T g is equivalent, in the sense of "universality 
classes" , to that of Ising spin glasses in a magnetic field 
h @. Spin-glass droplet scaling ideas [!, 0] were used 
to discuss the consequent expected glass phenomenology. 
This approach would be appropriate if the lengthscales 
of cooperatively rearranging regions near T g were many 
times larger than the intermolecular separation. In fact, 
recent studies Q have indicated that this lengthscale is 
rather modest and only a few intermolecular distances. 
As a consequence, glasses are not really in the regime 
where droplet scaling ideas are appropriate. We shall 
therefore examine in this paper the glass phenomenology 
which arises when the correlation length is not large but 
instead is in the precursor regime to the droplet scaling 
limit. Rather to our surprise, we are able to find in this 
regime all the characteristic features of glass phenomenol- 
ogy such as the Vogel-Fulcher relation [l[ , the Kauzmann 
temperature 0], the Adam-Gibbs relation 0, etc... 

We shall study in particular the Edwards-Anderson 
(EA) Ising spin glass model @ in the presence of an ex- 
ternal magnetic field h both in one dimension (Id) and 
in three dimensions (3c?). The behavior in both dimen- 
sions is similar, but the Id case can be studied more 
thoroughly as its equilibrium properties can be deter- 
mined exactly by means of a renormalization group ap- 
proach, and its dynamical properties are accessible via 
Monte Carlo simulations. Even in the Id case, the model 



is able to mimic most of the experimental observations 
on supercooled liquids. An apparent Kauzmann para- 
dox [(J is found, accompanied by a growing (but still 
modest at the fields which we use) lengthscale £ and 
by an apparent divergence of the relaxation time as in 
the Vogel-Fulcher (VF) equation with T VF = T K 
There is thus an apparent thermodynamical and dynam- 
ical singularity at a finite temperature Tk , but it is not a 
true transition. Tk is just a crossover temperature such 
that when T < Tk the growth of £ as T decreases has 
largely ceased. The 3d case cannot be solved exactly but 
has been studied within the Migdal-Kadanoff approxima- 
tion (MKA) [l(| , and a similar glass phenomenology 
emerges. But there are some significant difference be- 
tween the Id and 3d cases, due to the fact in Id there is 
no finite temperature spin glass transition, but only a di- 
verging lengthscale as T — > when h — 0, whereas in 3d, 
there is in zero field a finite temperature transition and 
so in order to have a lengthscale of only a few intermolec- 
ular diameters at low temperatures, a large field has to 
be applied. Our work within the MKA in 3d does not 
provide a quantitatively accurate picture, but the results 
are so encouraging that it would seem worthwhile to at- 
tempt to get more quantitative results, probably by use 
of Monte Carlo methods. Unfortunately glass timescales 
are so long compared to molecular collision timescales 
that realistic simulations will be challenging. 

The investigations presented here suggest a new frame- 
work to understand glass behavior: features of the inter- 
molecular potential and the density determine the value 
of the field h and the temperature scale, but once these 
are fixed there is sufficient universality left near the glass 
transition T g that the mapping to Ising spin glasses in 
a field provides a semi-quantitative account of both the 
thermodynamic properties of glasses and those dynami- 
cal features which can be undestood in terms of flipping 
and cooperative rearranging of spin domains of linear ex- 
tent £ sitting in a random effective magnetic field, requir- 
ing free energy activation over barriers. 

This paper is organized as follows: in Sec. [IT] we solve 
the EA Ising spin glass in an external magnetic field 
in one dimension by using a decimation approach and 
we study its dynamical properties by performing Monte 
Carlo simulations; we then discuss the connection to the 
phenomenology of glasses. In Sec. Mil we examine the 
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three dimensional case, by means of the Migdal-KadanofT 
approximation. Finally, in Sec. HVl we present some con- 
cluding remarks. 



II. THE EA ISING SPIN GLASS IN ONE 
DIMENSION 

The EA Ising spin glass Hamiltonian in Id in the pres- 
ence of a uniform external magnetic field reads: 



h = - cr.i+1 - k o % , 



(i) 



where the spins <jj can take values ±1, and the nearest- 
neighbor couplings Ji are independent of each other and 
Gaussianly distributed with zero mean and standard de- 
viation J. In principle h could be a function of tempera- 
ture but we shall regard it as a temperature independent 
constant, whose magnitude is chosen both in lei and 3<i 
so that the low-temperature spin glass correlation length 
£ is of the order of a few lattice spacings and so is com- 
parable to the glass correlation length of real glasses at 
T g We evaluate the free energy of the system by 
using an iterative real-space renormalization group tech- 
nique (9l.[lo|. It consists in tracing out every other spin in 
the system, thereby generating new effective interactions 
between the remaining spins which sit in new magnetic 
fields: 

ST e P[4 n) ^+i+4ll^+i^+2+J2]^ h( r' > ^] ( 2 ) 



At the n-th step in the decimation process, 4 n) and h\ n) , 
have probability distributions, which evolve with the it- 
eration. In Id, the EA model has a genuine critical point 
at (T = 0, h = 0), corresponding to a non-trivial fixed 
point of the recursion relations. Conversely, at any finite 
temperature and magnetic field, the system evolves to- 
ward a trivial "random paramagnetic" fixed point: the 
variance of the effective couplings decreases under itera- 
tion and approaches zero, whereas the effective magnetic 
fields have a distribution which approaches a Gaussian, 
with mean h (i.e., the initial value of the magnetic field), 
and variance <Jh(T, h). 

For each realization of the quenched disorder, the free 
energy, fj, can be determined exactly by summing the 

(n) 

spin-independent terms, W i __, which are generated at 
each step of the decimation Once the average over 
the disorder is taken, / = [fj]j, the entropy density is 
obtained using S = —df/dT\h- In Fig.QJi the entropy is 
plotted versus the temperature for three different value 
of the magnetic field, h = 0.05, 0.125 and 0.2. The fig- 
ure shows a temperature range in which the entropy de- 
creases linearly, and would be extrapolated to vanish at 
T K (h) as 



demonstrating that the model has a "Kauzmann para- 
dox" similar to that observed in supercooled liquids. 
However, below a crossover temperature, T*(< T g ), the 
entropy deviates from linearity and does not vanish com- 
pletely except at T = 0. 

By computing the derivative of the free energy with 
respect to the variation of the magnetic field on two dif- 
ferent sites, it is possible to evaluated the equilibrium 
connected correlation function: 



(a i( T i+l } 2 c = T 2 



d 2 InZ 



d{5hi)d{5h t+l ) 



(4) 



8hi=8h i+l =0 / 



^From the exponential decay of the correlation function, 
[(iTj(Ti+z)c] j ~ exp(— Z/£), we extract the equilibrium cor- 
relation length of the system, £, plotted in Fig.QJ) for the 
same values of the magnetic field as before. £ increases as 
the temperature is decreased but at low enough temper- 
atures (T < T*), it bends over and approaches a finite 
value at T = 0, proportional to /i -2 / 3 , as predicted by 
the droplet picture [3, [J] on equating the energy to flip a 
domain of size £ to the magnetic field energy which could 
be gained, J£ e - h£, d ' 2 {~ k B T K ), and for d = 1, 9 = -1. 
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S~k B c{h)[T-T K {h)], 



(3) 



FIG. 1: One dimensional EA model in a field for h — 0.05 
(continuous line and circles), h — 0.125 (dashed line and di- 
amonds) and h — 0.2 (dashed-dotted line and triangles), a) 
Entropy per spin, S; b) Equilibrium correlation length, £; 
c) Adam-Gibbs relation: temperature dependence of S£; d) 
Logarithm of the relaxation time, t, as a function of the tem- 
perature. The curves correspond to the Vogel-Fulcher fits, 
Eq. ©, with Tvf = Tk- The arrows indicate the Kauzmann 
temperatures, Tk, whereas the vertical dashed lines corre- 
spond to the crossover temperatures, T* . 

In order to establish a connection with the dynamical 
features of glass-forming liquids, we have studied the dy- 
namical properties of the Id model by performing Monte 
Carlo simulations of a system of 1024 spins. We have 
computed the spin-spin cumulant auto-correlation func- 
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tion, defined as: 



C(t, t w ) 



(5) 



For large enough waiting times, t w , the system reach sta- 
tionarity, characterised by time transitional invariance, 
i.e., C(t, t w ) — C(t). Although the mapping established 
in Ref. [2J is explicit only for equilibrium quantities and 
might not extend to dynamical features, we find that, in 
analogy with glass-formers, C(t) is very well fitted by a 
stretched exponential form, C(t) ~ exp[— (t/r) 13 ], where 
r(T, h) is the system relaxation time. The spin-spin auto- 
correlation function is plotted in Fig. [2] for several values 
of the temperature. 

The relaxation time is plotted in Fig. [TJi for h — 0.05, 
0.125 and 0.2 as a function of the temperature. Simi- 
larly to what happens in supercooled liquids pj a Vogel- 
Fulcher law, 



log r = log T 



VF 



T — Tvf ' 



(6) 



is able to fit quite accurately the data for each value of 
the magnetic field (over 4-5 decades), with the VF tem- 
perature, Tvf, set equal to the Kauzmann temperature, 
Tk [H- Again, the dynamical singularity is only appar- 
ent, since the relaxation time diverges only at T = 
according to an Arrhenius law. At low enough tem- 
perature r deviates from the VF law. Such a depar- 
ture starts to emerge for h = 0.2 when T < 0.25. In- 
terestingly, the onset of the deviation from the VF fit 
seems to coincide with the crossover temperature T* at 
which the entropy deviates from linearity. This behav- 
ior is consistent with the Adam-Gibbs (AG) relation Q, 
logr = A ag + B AG /TS(T), which holds in the tem- 
perature range explored by the simulations, as shown in 
the inset of Fig. [5] The original derivation Q of the 
AG equation relies on the assumption of the existence 
of correlated regions of size £ rearranging cooperatively, 
and on the hypothesis that S£ d ~ const. As shown 
in Fig. [TJ;, this relation seems to hold in the interme- 
diate temperature window explored in the simulations, 
whereas it breaks down below T*. 

One might wonder whether the whole entropy, 5, 
which we have calculated, plays the role of the "configu- 
rational entropy" for supercooled liquids in the AG rela- 
tion. When T < Tk, the thermodynamics is dominated 
by the flipping of the few single spins for which the local 
field is comparable to T. For T > Tk, the thermodynam- 
ics is dominated by the excitation of "droplets" of size £. 
Hence it is tempting to regard the very low temperature 
entropy in our spin glass simulations as being the contri- 
bution to the entropy from a single state, and that the 
configurational entropy is that which arises when many 
droplets are thermally excited. We have studied also the 
"configurational entropy" which is obtained by subtract- 
ing from S the linear contribution which fits the very 
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FIG. 2: Main frame: Time dependent spin-spin auto- 
correlation function, C(t), ELS £L cL function of log t for h = 0.05 
and T = 0.45 (squares), T = 0.35 (triangles), T = 0.275 
(circles) and T — 0.225 (diamonds). The data are averaged 
over 16-32 independent realisations of the disorder and over 
t w . The dashed line correspond to stretched exponential fits 
C(t) ~ exp[-(t/r)' 3 ] with 0.2 < /3 < 0.3. Inset: Adam-Gibbs 
relation: Logarithm of the relaxation time, r, as a function 
of 1/TS for h = 0.05 (circles), h = 0.125 (diamonds) and 
h = 0.2 (triangles). The straight lines are guides for the eye. 

low temperature entropy. However, the same qualitative 
results were found with this definition of S. 

Due to the success of the AG relation, it is natural 
to expect that the dynamics of the system will be domi- 
nated by the flipping of spin domains of size £. Since the 
variance of the the effective couplings decreases under it- 
eration and approaches zero after a few iteration steps, 
whereas the effective magnetic fields are Gaussianly dis- 
tributed with mean h and variance ah (T, h) , one might 
guess that the dynamics of the system is equivalent to 
that of a chain of non-interacting spin domains of lin- 
ear extent proportional to £, sitting in a random exter- 
nal field. Notice that this situation has been studied in 



Ref. 11[ . Taking into account the time to pass a spin flip 
through the domain, which involves the breaking of the 
largest bond in the domain whose magnitude will be de- 
noted by Li, the time to reverse each spin domain will be 
of the form r oc exp[(2L^ + 2hi)/T}. If the distribution 
of the Li is also Gaussian, the distribution of the sum 
(Li + hi) will be another Gaussian of variance a 2 L + cr^T. 
According to Ref. [ll[ , this leads to the following expres- 
sions for the relaxation time r 
for the stretching exponent (3 ci 



-exp[4(of i + al)/T 2 } and 
C[l+4(al+al)/T 2 }-^, 
with a constant C = 1. We have verified that these for- 
mula work quite well in describing the dynamics of the 
system in one dimension (with <tl ~ 0.22). Nevertheless, 
they are not perfect: the constant C is bigger than one 
(C ~ 1.5) and <jl is too small. These discrepancies could 
be because the time it takes to flip the spins by breaking 
the largest bond on a line of spins of length £ has not 
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FIG. 3: Main frame: Kauzmann temperature Tk (continu- 
ous line and circles), jump in the specific heat AC P (dashed 
line and squares), and stretching exponent f3 (dotted line and 
triangles) as a function of the fragility D extrected from the 
VF fit of logr in Id. Inset: Kauzmann temperature Tk 
(continuous line and circles) and jump in the specific heat 
AC P (dashed line and squares) as a function of the fragility 
D extracted from the VF fit of exp[4o-^/T 2 ] in 3d. 

been handled with sufficient accuracy. The largest bond 
has its own probability distribution, which is just not a 
Gaussian. A full theory would be complicated. However, 
in 3d much bigger values of the external field have to 
be taken to keep the magnitude of the correlation length 
only a few lattice spacings at low temperatures, so one 
could reasonably expect that the relaxation time is dom- 
inated just by the random fields alone. 



0.3 






0.2 






0.1 




a 













3 




b 


2.5 












2 






1.5 







12 3 4 5 

T 



FIG. 4: Three dimensional EA model within the MK approx- 
imation for h = 2.0. a) Entropy per site, 5"; b) Correlation 
length, £. The vertical arrows and the vertical dotted lines 
correspond, respectively, to Tk and T*. 

The fragility D, and the Kauzmann temperature Tk 



(see Eq. (|6])) are both affected by the magnetic field, 
h. More precisely, D decreases as h is increased (i.e., 
stronger glasses are described by bigger fields), whereas 
Tk (and T*) increases as h is increased. In the main 
frame of Fig. [3l the dependence of Tk on the fragility is 
plotted. We also plot the behavior of another two impor- 
tant quantities: the jump in the specific heat AC P and 
the stretching exponent j3 (at T = 0.3). The former can 
be estimated by (arbitrarily) setting T g ~ T*, so that 
AC p ~ c(h)T*, c(h) being the slope of the entropy in 
the linear regime [see Eq. ([3])] ; the latter can be directly 
computed from numerical simulations. 

Interestingly enough, f3 is found to be an increasing 
function of D, whereas AC P decreases as D increases, in 
agreement with observations on supercooled liquids [12j . 



III. RESULTS IN THREE DIMENSIONS 

We now turn to the 3d case, which is most relevant 
for real supercooled liquids. We have evaluated the free 
energy of the EA spin glass model in a field by means 
of the MK approximation, a real-space renormalization 
group techniques that gives approximate recursion rela- 
tions for the flow of the coupling constants and magnetic 
fields distributions @, [lj| • We have used the "bond mov- 
ing" procedure, where the bonds on the 3c? lattice are 
moved before each decimation step, so that no higher or- 
der couplings are generated [io| : in a d-dimensional lat- 
tice, 2 d_1 bonds are superimposed and added up, whereas 
the "naked" spins that are left behind have no cou- 
plings. Taking the trace over the spins that are on the 
main bonds leads to the coupling constants, according to 
Eq. $2$, between neighboring spins on the coarse-grained 
lattice. The decimation procedure is iterated n times on 
a lattice of size L = 2™. There is a flexibility in the MK 
renormalization scheme as to how the fields are moved. 
We have treated the field terms as belonging to bonds: 
when a bond is moved we also move all its field terms to 
the end that is to be traced over next pL0j | . 

In 3d, within the MKA, when h = and T < T c , 
T c ~ 1.78, the variance of the effective couplings grows 
indefinitely under iteration. For finite values of the mag- 
netic field there is no evidence of a de Almeida-Thouless 
line [III]. The variance of the couplings might grow ini- 
tially for low enough fields and/or temperatures, but it 
always decreases and eventually vanishes after a sufficient 
number of iteration steps, just like in Id. The average 
value of the field distribution equals the initial value of 
the uniform magnetic field, h, whereas the width of the 
distribution saturates at a finite value, cr^(T, h). 

The temperature dependence of the entropy per spin 
and the correlation length £, obtained within the MKA 
from the exponential decay of the variance of the effec- 
tive coupling, which decreases as exp(— 2™/£) at 
large n, are plotted in Figs.Hk and 2b for h — 2.0, show- 
ing a scenario very similar to that found in lei: there is 
a temperature range in which the entropy per site, S, 
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is linear and is extrapolated to vanish at a finite Kauz- 
mann temperature, TV; a crossover occurs at a higher 
temperature T*, where the entropy deviates from linear 
behavior. The correlation length, £, grows as the tem- 
perature is decreased and approaches a finite value at 
T = 0. 

We also mention that, similarly to the Id case, there 
is a modest range of values of the external magnetic field 
h (for 1.8 < h < 2.5) for which the AG relation, ~ 
const, holds in the temperature region T > . This AG 
relation, however, breaks down below T* and at high 
temperatures and in contrast to the Id case is of less 
utility. 

Due to the magnitude of the relaxation times, standard 
Monte Carlo simulations of the 3d model are more chal- 
lenging than in the Id case, and we will leave them for 
future investigations. However, since the values of inter- 
est of the magnetic field are much bigger than in Id, one 
can surmise that the flipping of a spin domain of size £ sit- 
ting in a random external magnetic field Gaussianly dis- 
tributed and with variance a 2 are the dominant dynami- 
cal processes; hence, the energy barriers involved in such 
processes might account reasonably well for the system's 
relaxation time, leading to logr ~ 4er 2 /T 2 [11|. Follow- 
ing this hypothesis, we have verified that exp[4cr 2 /T 2 ] can 
be well fitted by a VF law with T VF = T K for T > T*. 
From the VF fit of this quantity it is also possible to ex- 
tract the fragility, D, in the 3d case. In the inset of Fig.[3j 
Tk and AC P are plotted as a function of D, showing they 
have very similar behavior to that found in Id. 

The domain size £ according to droplet scaling is given 
by equating the cost of flipping a droplet of size £, £ e , 
to the field energy which might be gained, h£ d / 2 . As in 
3d the exponent 9 is small, (« 0.2) it follows that 
a h ~ h 2 £f ~ const., which implies that the AG relation 
logr ~ A/TS should hold. 

The stetching exponent would be expected to be [1 + 
4<t 2 /T 2 ]~ 1 / 2 , provided again that the time taken to pass 
the domain wall through the domain is not significant. 



One feature of the MKA study in 3d is that the 
configurational entropy seems to be smaller than the 
quoted values near the glass transition, perhaps by as 
much as a factor of 3 12]. In the derivation of the 
mapping to spin glasses [2| one can see that the field h 
will be a function of both the temperature and density, 
rather than simply being a temperature independent 
constant as we have assumed here throughout for 
simplicity. Allowing for this temperature dependence 
could significantly change the entropy. For example 
suppose h 2 — hg + b 2 T 2 , then the high-temperature limit 
of the entropy is S = ln[2cosh&], and by adjusting b, can 
be made as large as desired. 

IV. CONCLUSIONS 

The mapping between supercooled liquids and spin 
glasses in an external magnetic field, proposed in Ref. [2], 
thus seems to provide a semi-quantitative explanation of 
the properties of supercooled liquids including the Kauz- 
mann paradox, the Vogel-Fulcher behavior of the relax- 
ation time, the stretched exponential decays of corre- 
lation functions, a growing lengthscale, and the Adam- 
Gibbs relation in the regime T > T g , which is the precur- 
sor regime accessed by the experiments, where the cor- 
relation length is growing with temperature but is still 
only a few intermolecular distances. The droplet scal- 
ing limit studied in Ref. [2| is appropriate only when the 
correlation length is much bigger than this. The large 
timescales which exist below the glass transition temper- 
ature T g prohibit the taking of equilibrium data below 
it and so the apparent thermodynamic and dynamical 
singularities at Tk cannot be accessed. In our scenario, 
Tk is only a crossover temperature at which the growing 
correlation length saturates to a constant value. 
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